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Abstract. Shifts from the expected nuclear magnetic resonance frequencies of 
antimony and bismuth donors in silicon of greater than a megahertz are observed in 
electrically detected magnetic resonance spectra. Defects created by ion implantation 
of the donors are discussed as the source of effective electric field gradients generating 
these shifts via quadrupole interaction with the nuclear spins. The experimental results 
are modeled quantitatively by molecular orbital theory for a coupled pair consisting of 
a donor and a spin-dependent recombination readout center. 
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Group-V donors are the most studied impurities in silicon (Si). For decades, the 
shallow donors have been employed for the fabrication of transistors in microelectronics. 
More recently, they have become very promising candidates as quantum bits [1], since 
they can be considered as hydrogen-like atoms embedded in a solid-state matrix. The 
weak spin-orbit coupling present in the silicon host makes the coherence time of the 
donor electron and nuclear spins exceptionally long [2, 3]. Moreover, the existence of 
spin-free 28 Si allows for the isotopic purification of Si [4, 5, 6], further improving these 
coherence times by removing most of the 29 Si (nuclear spin / = 1/2) [7, 8, 9, 10, 11]. 

This has motivated researchers to investigate techniques to incorporate donors in silicon 
to produce single donor devices [12, 13, 14]. The electron and nuclear spin (/-factors of 
shallow donors are already well-known [15]. Similarly, the hyperhne interaction between 
the donor electron with donor nuclear spin [15, 16] and with 29 Si nuclear spins [17] have 
been measured and are used to understand the donor electron wavefunction [18, 19]. 
More recently, the Stark shift of the donor states has been observed [20, 21, 22] and 
modeled [23, 24, 25, 26, 27], 

Beyond their coherence properties, group-V donors in silicon with large nuclear 
spins (> 1/2) such as 75 As, 121 Sb, 123 Sb, and 209 Bi (see Table 1) provide 

additional accessible states which are very attractive resources for performing quantum 
computations. Muthukrishnan and Stroud [28] have shown that the use of qudits [29], 
d-level systems, allows to decrease the number of necessary elementary qubits for a 
computation by a factor of log 2 d. Moreover, Lanyon et al. [30] have demonstrated that 
the five two-qubit gates required to realize a Toffoli gate (a three-qubit entangling gate) 
can be reduced to only three gates if at least one of the three qubits has an additional 
available level. Also, an implementation of the Grover algorithm directly in large nuclear 
spins in semiconductors has been proposed by Leuenberger and Loss [31]. 

This paper reports the experimental observation of line shifts of the nuclear spin 
resonance of 121 Sb, 123 Sb and 209 Bi donors in Si at low magnetic fields observed via the 
electrical detection of mixed spin states. These shifts are interpreted as resulting from 
the nuclear quadrupole interaction (abbreviated NQI) of the nuclear electric quadrupole 
moment with the donor electron wavefunction, which is modified by the presence of 
implantation defects. A shift of the nuclear magnetic resonance (NMR) due to NQI 
arises from the interaction of a nuclear spin I > 1 with an electric field gradient (EFG) 
and is well-studied in the atomic physics and nuclear magnetic resonance communities 
[32, 33, 34]. However, the observation of a quadrupole-induced magnetic resonance shift 
of the group-V donors in Si, where the NQI with crystal field gradients vanishes due 
to the cubic symmetry, has only recently been realized for ionized [35] and neutral [36] 
donors in strained Si, and for 209 Bi implanted silicon devices [37]. While an external 
application of EFG large enough to result in an observable shift (> 10 18 V/m 2 ) is 
challenging, gradients of such a magnitude can be generated by donor electrons whose 
spatial symmetry is broken by the presence of strain [38], close acceptors [39] or ion 
implantation defects [40, 41, 42], As our detection method is based on the spin- 
dependent recombination (SDR) [43, 44, 45] of close pairs (1 ~ 10 nm) formed by 
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the donor electron and such defects, only donors that are exposed to these EFG are 
probed. By working at low magnetic fields, where the donor electron and nuclear spin 
wavefunctions mix strongly due to the hyperhne interaction, we are able to use such 
SDR signal to investigate possible quadrupole-induced NMR shifts [44], 

Three samples, Si: 121 Sb, Si: 123 Sb, and Si: 209 Bi, were employed. The substrate used 
for fabrication was from a highly resistive (> 3 kfhcm) float-zone Si wafer. The samples 
were ion-implanted with 121 Sb, 123 Sb or 209 Bi at room temperature. Bi was implanted at 
300 and 550 keV with doses of 0.7 x 10 13 and 1.3 x 10 13 cm" 2 , respectively. The isotope- 
selective Sb implantations were performed with an energy of 30 keV and a fluence of 

1.7 x 10 11 cm" 2 . These conditions yielded a maximum Bi plateau-like concentration of 

1.8 x 10 18 cm" 3 [above the solubility limit [46]] in a depth of 90 to 150 nm from the 
surface, and a 121,123 Sb peak concentration of 1.1 x 10 1 ' cm" 3 24 nm from the surface. 

The three samples were individually annealed after the ion implantation, at 650 °C for 
30 min in an evacuated quartz tube. This process, designed to maximize the number 
of donor-readout center (D-R) pairs [47], led to an activation efficiency [48, 49, 50, 51] 

< 60%, resulting in a Bi donor concentration less than 1.1 x 10 18 cm” 3 , and a 121 Sb and 
a Si: 123 Sb concentration less than 4 x 10 16 cm" 3 (in all cases below the metal-insulator 
transitions [52, 53]). 

In the SDR detection technique, the photocarriers generated by a 100-W halogen 
lamp are captured by the ionized donors. The recombination time of the donor electrons 
depends on the relative electron spins of the donor and the paramagnetic readout 
center: in this D-R pair, the antiparallel electron spin configuration leads to a fast 
recombination, whereas for the parallel spins, the expected recombination time is much 
longer (spin blockade) [54]. In a steady state and for a continuous illumination, mostly 
the parallel electron spin pairs remain without external induction of magnetic resonance. 
Therefore, flipping the donor electron spins by magnetic resonance enhances the electron 
recombination, leading to a decrease in the sample photoconductivity from the steady 
state. This change in photoconductivity was measured by probing the change in the 
microwave reflectivity of the sample in a JEOL JES-RE3X X-band cw-electron spin 
resonance (ESR) spectrometer tuned at 9.086 GHz and 100 mW, with an Oxford 3 He- 
cryostat and a lock-in amplifier using a modulation amplitude of 7 G at 100 kHz. For the 
low frequency excitations, a dedicated coil was used. A static magnetic held is applied 
in the crystallographic [001] direction. 

The donor spin Hamiltonian considered is 

Ft = g e fi B B z S z - g n [i N B z i z + hAS ■ I + I ■ Q ■ I, (1) 

where g e and g n are the donor electron and nuclear ^-factors, respectively, A is the 
strength of the isotropic hyperhne interaction (Table 1), and B z is the magnetic 
held in the z direction. The eigenstates |i) of Eq. (1) are cos</>; \ms = 1/2,raj) + 
sin0j \mg — 1/2, mj + 1), where mj and ms are the eigenvalues of I z and S z , and 
is the spin mixing angle, explicitly described in Ref. [56]. They are labeled with 
integers 1 < i < (2 S + 1)(2J + 1) in order of increasing energy. The term I ■ Q ■ I 
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represents the quadrupole interaction of the donor nuclear spin. For an electrostatic 
field E z experienced by the nucleus in presence of the electronic charge configuration of 
its atom, the effective EFG in z direction is dE z /dz = V z e f, and the NQ1 is written as 
[32] 

1-Q-l =hv Q (3 F - 1(1 + 1)) /(4/(2J - 1)), (2) 

where the coefficient hu Q = eQV z e z s , e is the elementary electrical charge and Q is the 
quadrupole moment of the nucleus. 

At high magnetic fields, the eigenstates of Eq. (1) are given by the products of the 
electron and nuclear spin states. In this case, their magnetic resonances are addressed 
separately, either in ESR (A mj = 0, A ms = ±1) or NMR (A ms = 0, Am/ = ±1). 



Figure 1 . (Color online.) (a) ESR and (b) NMR transition frequencies and (c) NMR 
resonant magnetic field sensitivity to the NQI of Bi donors as a function of magnetic 
field. The colored lines (| 19) «->■ 120) in blue, 110) -<->■ 111) in green and |1) e->- |2) 
in red) have high SDR intensities and were analysed in this study. The same color 
code is used in Fig. 2 (a, b and d). (d) Experimentally observed SDR spectrum of 
the Bi donor spin transitions under 200-MHz excitation. The black lines represent 
the expected resonance position of the electron and nuclear spin transitions without 
exchange interaction or NQI. The arrows indicate the assignment of the spin transitions 
between the states | i) and | j), for an angle d = 0 or 0 = tt/2 (see text). 
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Table 1. Magnetic resonance parameters of group-V 75 As, 121 Sb, 123 Sb and 209 Bi 
donors in Si. The quadrupole moment of the donor nuclei Q is given in 10 28 m -2 . cn 
is the Bohr radius as defined in [27]. 



Si: 75 As 

Si: 121 Sb 

Si: 123 Sb 

Si: 209 Bi 

I 

3/2 

5/2 

7/2 

9/2 

9e 

1.99837 a 

1.9986 a 

2.00049(5) b 

9n 

0.9599® 

1.3454® 

0.7285® 

0.9135® 

A (MHz) 

198.35 a 

186.80 a 

101.52 a 

1475.05(17) b 

di (nm) 

1.45 d 

1.67 d 

0.967 d 

Q(b) 

+0.307(4)® 

-0.41(4)® 

-0.49(5)® 

-0.58(9)® 


a Feher (Ref. [15]) 
b Mortemousque et al. (Ref. [45]) 
c Wolfowicz et al. (Ref. [10]) 
d Pica et al. (Ref. [27]) 
e Stone (Ref. [55]) 


In SDR spectroscopy, which is sensitive to the electron spin configuration, only terms 
involving the electron spin operator S can be detected, namely the electron Zeeman 
interaction and the hyperfine interaction. For hA > g e nsB z , however, the eigenstates 
are linear combinations of several nuclear and electron states, so-called mixed states 
[57, 58]. Resonant transitions which change both spin states become possible, allowing 
for transitions which are sensitive to the nuclear spin properties but also change the 
electron spin state to become detectable by SDR [44, 59]. The sensitivity of the 
low-field line positions (resonant magnetic fields) to the NQI, noted dB§/dvQ, were 
calculated numerically for 209 Bi [Fig. 1(c)]. The spin-dependent conductivity signal is 
based on coupled spin pairs, which are most efficiently formed for ‘pure’ electron spin 
states [57, 44], The transitions most easily detected [Fig. 1 (b, c)] are therefore those 
involving the two product states: \ms = +S) \mj = + 1) = 120) (119) -B- 120), in blue) 
and | ms = —S) \mi = — I) = 110) (110) -H- 111), in green). Other NMR transitions have 
a weaker SDR intensity [Fig. 1(d)]. 

The ESR and NMR transition frequencies of Bi donors are plotted in Fig. l(a, 
b). Figure 1(d) shows an SDR spectrum of 209 Bi donors in Si recorded under a 200- 
MHz excitation. Four of the observed peaks are identified as the resonance transitions 
induced by the X-band microwave employed for the SDR detection of the ac-conductivity 
of the sample. They can be used to extract the electron (^-factor and the hyperfine 
interaction A since they are, in good approximation, not influenced by the nuclear 
Zeeman interaction or the NQI. As discussed in Ref. [45], the observed values are 
slightly different from those observed in conventional ESR due to the influence of the 
recombination partner. In addition, several peaks are observed due to the rf excitation, 
in particular two prominent lines at ~ 575 G and ~ 960 G attributed to the NMR 
transitions 110) -H- 111) and 119) -H- 120), respectively. However, when comparing these 
to the line positions as calculated for Eq. (1), using the experimental values of A and g e 
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extracted from the X-band line positions and assuming zero NQI [Fig. 1(d)], the NMR 
lines appear to be shifted significantly. These line shifts could not be observed in our 
previous report [44] due to the lower concentration in Bi donors. We could not observe 
forbidden electron spin transitions (except |1) yy 1 19)) or spin transitions arising from 
mixed mw-RF excitation frequencies. Octupole interaction was not observed and higher 
order multipole interactions are assumed to be negligible. 

To gain further insight into the origin of the observed shift, the resonance lines are 
recorded at different excitation frequencies as shown in Fig. 2(c), where the horizontal 
axes are aligned to the expected resonance held B Q h for the spin system without NQI. 

For all the measured NMR spectra of Si: 209 Bi, the shift of the resonant magnetic held 
for three transitions is shown in Fig. 2(a) as a function of the magnetic held. Since the 
observed shifts are orders of magnitude larger than could be expected for a change in 
nuclear or electron g-factor, or in hyperhne interaction, we will now discuss the NQI as 
a possible explanation. To separate the effect of the NQI from other possible influences, 
the data are plotted against dB 0 /duQ of the corresponding transitions [Fig. 2(b)]. The 
line positions extracted from the 1 19) yy 1 20), 1 10) yy 1 11) and |1) yy |2) transitions are 
well described by a linear dependence, which confirms the NQI of the nucleus as origin 
of the observed shift. The linear ht slopes of Fig. 2(b) indicates a NQI v Q e fff = —18.2(2) 

MHz in Si: 209 Bi. The corresponding effective EFG is = 1.30 x 10 21 V/m 2 . 

However, an additional offset is necessary to reproduce the observed line shifts. 

A possible explanation is the exchange interaction {hJ S R ■ S ) between the readout 
partner (electron spin state m R ) and the donor electron spins which was not considered 
in Eq. (1), and leads to a constant shift of the resonance for the transitions and 
frequency region considered here. This is supported by the different offset signs for 
the three transitions, in agreement with the steady state polarization of the spin 
pairs due to the spin-dependent recombination process. With the expected sensitivity 
dBo/dJ = -m R j^p- & —0.18 G/MHz for 1 19) yy |20) (+0.18 for 1 10) yy 1 11) 
and |1) yy |2)), the observed +8.0(4) G offsets in Fig. 2(b) can be explained by a 
coupling constant J(Bi) = 45(2) MHz, significantly stronger than usually observed in 
SDR experiments [60]. Another effect contributing to these offsets is the asymmetric 
and broad distribution of the NQI (assumed exponential). Such a distribution 
makes the resonance line asymmetric (see Fig. 2(c)). The line shape, symmetric 
for dBo/du Q ps 0, becomes very asymmetric at more negative dBo/du Q , revealing 
the NQI distribution. More quantitatively, the parameter 7+ y of the fitting function 
d 2 {exp(— B/rasy) G> exp(— (B — B exp ) 2 /2cr 2 )} /dB 2 is plotted in Fig. 2(d). The line 
offsets in the latter indicate that the asymmetry of the resonance peaks is not only 
caused by NQI, but by exchange interaction as well. 

To investigate further the origin of the EFG responsible for the NQI, we reproduced 
the same measurement protocol with the two isotopes of the Sb donors in Si, which have 
an equal Bohr radius, different from Bi donors, allowing us to characterize the role of 
the donor electron distribution in the generation of the EFG. The shift of the resonant 
magnetic held for the transitions of the two samples are shown in Fig. 3(a) as a function 
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B z -B%( G) dBoldVQ (G/MHz) 


Figure 2. (Color online.) Resonant magnetic field shift Bg Xp — B$f of 119) -O- 120) (blue 
crosses), j 10) -<->■ 111) (green crosses), and |1) <->■ |2) (red crosses) of Si: 209 Bi plotted 
as a function of B z (a) and the sensitivity dB 0 /dv Q (b). The circles represent the 
resonant magnetic field shifts estimated from Fig. 1(d), but not taken into account for 
the fit of the NQI and exchange interaction parameters. The black crosses correspond 
to the |1) «->■ 120), 12) «-»■ 119) and |3) «->■ 118) X-band ESR transitions. Bg Xp ' and B}f 
are the measured and calculated resonant fields (assuming that v Q =0). The dashed 
lines are calculated from a single linear fit (see text), resulting in v Q = —18.2(2) MHz 
and an exchange interaction J = 45(2) MHz. (c) Recorded SDR spectra of the Bi spin 
transitions 119) <->■ 120) (blue) measured at various frequencies. The spectrum (black) of 
the 19) -<->■ 112) transition at the hyperfine clock transition [45] is plotted for comparison 
(insensitive to distribution in hyperfine interaction, and very weakly sensitive to NQI). 
(d) Asymmetry parameter r asy (dB 0 /du Q ) (see text) plotted for Si:Bi spectra (with 
high enough signal to noise ratio and non overlapping transition lines). The dashed 
lines (a and d) are guides for this eye and were calculated assuming u Q = —18.2 MHz 
and J = 45 MHz. 


of the magnetic field. Figure 3(b) [respectively (c)] shows the measured shift of NMR 
transitions of 121 Sb [ 123 Sb] donors against dB 0 /dv q . A linear fit of the points indicates 
a NQI z/® xp ( 121 Sb) = —2.90(7) MHz and iA xp ( 123 Sb) = —2.95(9) MHz, and an exchange 
interaction of J( 121 Sb) = 13(2) MHz and J( 123 Sb) = 20(3) MHz. 
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The quadrupole moment ratio Q^isb/Q^ssb is 0.85(17). An equivalent ratio is 
expected for i/® xp ( 121 Sb)/z/® xp ( 123 Sb). In this study, the experimental value of this ratio 
is 0.98(5), which is in good agreement with the theoretical value. Though further 
studies could help to quantitatively explain the measured line positions, e.g. with 
electron nuclear double resonance experiments in the high-field limit, our experiments 
indicate NQI in the order of several MHz, corresponding to effective EFG on the order 
of 10 2 ° V/m 2 . 

The NQI with an uniaxial EFG shows a dependence on the angle 6 between 
gradient axis and the quantization axis [61], which is in first order (3 cos 2 6 — l)/2. The 
tetrahedral symmetry of the donor is broken by the presence of a defect center and for 
a random distribution of defect sites, one expects a continuous distribution of eQV zz /h 
from —VQ ax /2 to Uq ax for a quantization axis independent from the NQI. The lock-in 
detection of the change in the sample conductivity used in this study allows us to detect 
only the edge of the broadened resonance peaks. Therefore a continuous distribution in 
V zz leads to an effective peak splitting [blue labels in Fig. 1(d)]. The signal-to-noise ratio 
achieved experimentally did not allow the observation of all resonance peak splittings 
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Figure 3. (Color online.) (a) Detuning B g Xp — Bg h of the transitions 111) -O- 112) 
in Si: 121 Sb (blue) and 115) O 116) in Si: 123 Sb (red) plotted against B}f. The 
corresponding lines are guides for the eye calculated using the fitting parameters, (b, 
c) Shift, as defined in Fig. 2, observed for various transitions of Si: 121 Sb (b) and 
Si: 123 Sb (c) plotted as a function of dB 0 /dv Q . The linear fits indicate ^ Q ( 121 Sb) = 
-2.90(7) MHz and i/ Q ( 123 Sb) = -2.95(9) MHz. 
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Table 2. Simulated values, using Eq. (3), of the EFG for different valley configurations 
for I > 1 group-V donors. Only values of EFG larger than 10 15 V/m 2 are listed. 
The EFG numerical values are given in units of 10 19 V/m 2 . We use the valley 
notation Ai: (X + X + Y + Y + Z + Z)/\/6 , E x>y : (X + X — Y — Y)/ 2, and E 2z : 
(2Z + 2Z - X - X - Y - Y)/VV2. 





Si: As 

Si:Sb 

Si:Bi 

<Ai| 

'll orb 1 

'*'XX, | 

E x,y) 

3.0 

2.7 

11.6 

<Ai| 

nj orb 1 

7X yy. 1 

E x,y) 

-3.0 

-2.7 

-11.6 

(Ai| 

'17 orb 1 
^zz | 

E x,y) 

0.0 

0.0 

0.0 

(Ai 

1 ' 1 /orb 

| '‘•XX 

| E 2z) 

-1.7 

-1.6 

-6.7 

(Ai 

I nj orb 

1 rL yy 

| E 2z) 

-1.7 

-1.6 

-6.7 

(Ai 

1 'l/orb 

| '‘"ZZ 

j E 2z ) 

3.4 

3.2 

13.4 


and the spectral density of the resonances limited the peak assignment. 

No control of the electric field or of the EFG was performed during the experiments 
and the EFG induced by local strain [35] is negligible as no evidence of strain could 
be found in Si:Bi [45] nor in Si:Sb at X-band. Therefore, we propose that the EFG 
at the nucleus is induced by a redistribution of the donor electron density. To confirm 
this hypothesis, we proceeded in two steps. First, the orbital of the donor-readout 
center electron pair can be modeled to evaluate the mixing of the different valleys of 
the donor Is ground state. This was performed using a variational calculation using a 
wavefunction basis that includes the six donor valley configurations of the Is state plus 
a localized readout center [45]. The energy levels of the donor ground states are known 
from optical studies [62], whereas the localised readout center energy level is assumed to 
be E r = —0.25 eV [63]. The envelop functions were calculated using a corrected effective 
mass theory donor wavefunction. The parameter g = 159.4 described in Ref. [64] was 
used to adjust the envelop functions close to the donor nucleus [65]. Due to the presence 
of the readout center, the different valley populations of the donor ground state adjust 
to become if = ctAi0Ai + where 0* are the twofold E and threefold T 2 valley 

configurations. The coefficients cp depend strongly on the relative positions of the two 
recombination partners with respect to the crystal axes. 

The second step was the computation of the EFG for each donor valley 
configuration. Then, the EFG V XiXi generated by an orbital 0 was computed for each 
donor valley configuration as [32] 


<0|?C b J0> = 


4vre 0 e Si 


0 ' 


. 3x 7 — r z 


0 dr, 


( 3 ) 


where Xi = x, y , z. According to the group theory, the first order term of the EFG is 
of the form cp(A 1 |'H orb |0 i ), and is significant only for the E-doublet valley configuration 
(the second and higher order correction terms are negligible, see Table 2). 

For the SiBi sample, the average D-R separation has been calculated to be 2.3 di 
[45]. The estimated value of z/ h (Bi) for such D-R separation is 0.11 MHz, causing an 
EFG underestimated by a factor ^® xp (Bi)/zA h (Bi) = 160. This discrepancy might be 
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explained by the antishielding phenomenon. Upon application of an EFG, the energy 
levels of all shells of an atom are modified, giving rise to an EFG induced by the 
subshell electrons. Feiock and Johnson [ 66 ] have shown that a first order correction to 
the subshell energy levels gives the total EFG at the donor nucleus Vff 1 = (1 — 
where 7 is the antishielding or Sternheimer coefficient. The antishielding parameters 
of donors in Si have been measured only for limited species and ionization numbers. 

A donor forms four covalent bonds with the nearby Si atoms, and its core is therefore 
close to D 5+ plus the contribution of the four covalent electrons. For Bi donors [ 66 ], 
7 (Bi° + ) = —47.24 and the effect of the four covalent bonds would be to reduce this 
value down to 7 (Bi°) = —159. As Bi donor and Sb donors have different numbers of 
core electrons, the antishielding parameters may be very different. The determination 
of the antishielding parameters is beyond the scope of this work. Another hypothesis is 
the presence of several defects in the vicinity of a donor, leading to an underestimation 
of the orbital mixing. 

In summary, we have observed shifts greater than a megahertz of NMR transitions 
for Bi and Sb donors in Si, and their probable link to NQ1. The molecular orbital model 
used in this study points at the mixing of the donor orbital ground states Ai with 
the E doublet as a possible origin of the EFG. As a similar orbital state mixing can 
also be achieved with an electric held [27], or with a global [36] or a local [37] strain, 
the presented results will allow to relate, in the context of devices based on donors in 
silicon, the device physical properties (local strain, electric potentials defined by local 
gates) to the magnetic resonance properties of the donor spins. Finally, our results 
imply that by using local gates to modify and control the donor wavefunction, it would 
be possible to generate on-demand a local EFG at the donor nucleus which could be 
used to manipulate the donor nuclear spin. 
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